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(54) Tracking of sampling frequency in a DAB receiver 



(57) A DAB receiver in which, after a DBA baseband 
signal output from a receiving section is converted into 
a digital signal at a predetermined sampling frequency, 
the digital signal Is Fourier-transformed by a Fourier 
transformer, to be output has a sampling pulse genera- 
tor, a DFT window generator serving as a Fourier trans- 
formation execution timing, an impulse response 
arithmetic section for arithmetically operating impulse 
response of a transmission path on the basis of a Fou- 
rier transformation output of a phase reference symbol, 

FIG. 



an impulse barycenter calculator in the impulse 
response, and a sampling frequency control means for 
calculating a change amount of barycentric position on 
the basis of a position deviation between DFT windows 
in previous and current frames and barycentric posi- 
tions in the previous and current frames and controlling 
the frequency of a sampling pulse on the basis of the 
change amount of barycentric position. 
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Description 

BACKGROUND OF THE INVENTION 
5 1 . Field of the Invention 

The present invention relates to a receiver in digital audio broadcast (Digital Audio Broadcast: DAB) and, more par- 
ticularly, to a receiver in a digital audio broadcast in which 2N digital data are divided into N groups each having two 
bits, the first data of each group and the second data of each group are sequentially input to an Inverse Fourier trans- 

10 former as a real-number part and an imaginary-number part, respectively, the real-number part and imaginary-number 
part output from the inverse Fourier transformer are converted into analog signals, the analog signals are multiplied by 
cos and sin waves of a carrier frequency fc, respectively, the multiplication results are synthesized to be radiated in a 
space, signals radiated in the space are received, the received signals are multiplied by the cos and sin waves of the 
carrier frequency, respectively, the multiplication results are converted into digital signals at a predetermined sampling 

75 frequency, the digital signals are input to a Fourier transformer, and a real-number part arxl an imaginary-nunrt)er part 
output from the Fourier transformer are output as the first and second data. 

2. Description of the Related Art 

20 Digital audio broadcast (DAB) in which audio signals are converted into digital signals to obtain serial data, the 

serial data are divided into sets each consisting of 2N digital data, the 2N digital data are divided into N groups each 
having two bits, N carriers having different frequencies depending on combinations of 1 and 0 of each two bits are 4- 
phase-PSK-modulated, the modulated signals are frequency-multiplexed to be sent from a transmission station, the fre- 
quency-multiplexed phase-modulated signal is received by a receiver and demodulated to be an audio output is pro- 

25 posed. The DAB is being studied for practical use in Europe and the like. 

In this DAB scheme, in order to reduce influence of selectivity fading, pieces of information are divided in parallel, 
and modulation is performed by using a plurality of carriers (frequency interleave). As a result, even if any carrier 
receives fading, influence is reduced as a whole. Basically, the DAB scheme is a frequency division multiplex (FDM: 
Frequency Division Multiplex) scheme. In a simple FDM, an interval between carriers must be sufficiently large to avoid 

30 spectra from overlapping, and good frequency use efficiency cannot be obtained. For this reason, an OFDM (Orthogo- 
nal Frequency Division Multiplex) scheme is proposed. In this OFDM, carriers are arranged to satisfy orthogonal con- 
ditions and to allow spectra to overlap, good frequency use efficiency can be obtained, and IDFT (Inverse Discrete 
Fourier Transform) and DFT (Discrete Fourier Transform) operations can be used in a modulator and a demodulator, 
thereby advantageously obtaining very simple hardware. 

35 

(a) Principle of Digital Audio Broadcast by OFDM Scheme 

FIG. 13 is a view showing the theoretical arrangement of a transmitter for digital audio broadcast. Reference 
numeral 1 denotes a serial/parallel converter (S/P converter) for converting serial data d(n) (a(0), b(0). a(1), b(1 ),...) 

40 input at a transmission rate fs (= 2/M) into 2N-bit parallel data; 2q to 2n_i , N carrier multiplying sections for dividing the 
2xN bit parallel data into N groups a(0), b(0),; a(1), b(1);,..; a(N-1), b(N-1) each having two bits, multiplying thefirst bits 

a{0). a(1) a(N-1) of the groups by carriers (costont) having frequencies fo to fw-i. and multiplying the second bits b(0). 

b(1) b{N-1) by carriers (-sincont) having frequencies fo to fw-i *. and 3, a frequency multiplexer (MUX) for syrrthesizing 

output signals a(n)coso)nt and -b(n)sin<Ont (n = 0 to N-1) from the carrier multiplying sections of the respective groups 

45 and frequency-multiplexing the synthesized signal to send a signal D(t). 

When output signals from the carrier multiplying sections of the respective groups are synthesized by the frequency 
multiplexer 3, the carriers having frequencies fo to f|sj.i are 4iDhase-PSK-modulated by combinations of 1 and 0 of two 
bits of each group. The output D(t) from the frequency multiplexer 3 is given by: 

50 D(1) = X{a(n)cosci>nt - b(n)sinffl ^t} (n = 0 to N-1) 

Note that, when the frequency Interval of carriers Is represented by Af, the (n+1)th carrier frequency fn is given by: 

fn = f 0 + nAf 

55 

When a transmission lime of 2-bit data is represented by At (transmission rate fs = 2/ At), the frequency interval Af is 
given by: 
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Af = 1/NAt 



FIG. 14 is a view for explaining the function of the frequency multiplexer 3. N carriers fo to 1^^.^ having intervals Af 
are 44>-PSK-modulated by data a(0). b(0): a(1). b(1):...; a(N-1). b(N-1) of N 2-bit groups, and the modulated signals are 
frequency-multiplexed to be transmitted. FIG. 15 is a view for explaining a symbol. One symbol consists of 2xN bits. 
When the time length of one symbol is represented by Ts. the following equations can be obtained: 



Every symbol (2xN bits) Is subjected to 4<)>PSK modulation, and the modulated signals are frequency-multiplexed, 
thereby sequentially transmitting frequency multiplication signals D(t). 

FIG. 16 is a view showing a theoretical arangement of a receiver for digital audio broadcast. Reference numerals 
4o to 4^.1 denote N carrier multiplying sections for multiplying the received signal D(t) by carriers (cosont -sincont. n = 
0 to N-1) having frequencies fo to fisj.i; Sq to 5n-i. integrators for integrating outputs from the multiplying sections to 
demodulate data; and 6, a parallel/serial converter (P/S converter) for converting 2xN bit parallel data into serial data. 

The integrators 5o to 5|v|.i perform the following arithmetic operation to an input signal D(t): 



to demodulate data a(0), b(0); a{1), b{1);...; a(N-1), b(N-1). 

(b) OFDM Modulation Scheme using DFT 

In order to generate a baseband signal D(t) of OFDM, N 4(|>PSK modulators are required. In addition, N 4<|>PSK 
demodulators are required to demodulate the baseband signal D(t). When the number N is large, this scheme is not in 
practical use. Therefore, a method of simply performing modulation/demodulation by using DFT will be described 
below. 

(b-1) Anrangement of Modulator 

A baseband signal of OFDM is expressed by the following equation: 

[Equation 2] 



Ts = NAt 



Af = 1/Ts 



QEquation 1] 





N-1 




Jfc=0 



(1) 



When d(k) = a(k] -i- jb(k) is established, Equation (1) is expressed as follows: 
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[Equation 3] 

N-1 



-^RlJldikye'"'^] (2) 

Jfc=0 

70 where * means a complex number, and RQ represents the real -number part of Q. 
Since the following equation is established: 

[Equation 4] 
75 N-1 

N-1 

— Z [jiaj^+Jb^ cos ojj^t-h j(aj^+jbj^ sin coj^i] 

20 

N-I 

= (ajfoswj^t'bj^inujj^t) jXcj^ cj j^t bjfios cj f^l) 

25 =i?(O+y/(0 

it is apparent from Equation (3) that D(t) is expressed by Equation (2). In Equation (2), when I = mAt is established, 



[Equation 5] nli , ^ 



where it must receive attention that D(m) is the real-number part of ID FT (Inverse Discrete Fourier Transform) of d{k). 
When D(m) is output as 0(0), D{1),..., D(N-1) for each At, D(m) Is expressed as shown In FIG. 17. TTils signal must be 
obtained by sampling Equation (2) for each At. Therefore, when the signal passes through an ideal filter having the fre- 
40 quency characteristic shown in FIG. 1 8A and the inopulse characteristic shown in FIG. 1 86, the same equation as Equa- 
tion (2) can be obtained. 

As is apparent from the above description, in the modulator, ID FT is performed to input data d(k) (complex), and 
the real-number parts D(0) to D(N-1) are output for each At arxi are caused to pass through the ideal filter, thereby 
obtaining a modulation wave (baseband signal). FIG. 19 is a view showing the arrangement of a main part of a trans- 
45 mitter constituted while giving attention to this point. Reference numeral 1 1 denotes an IDFT section for performing 
inverse discrete Fourier transform to input data d(k) which is expressed as complex numbers by data a(k). b(k) (k = 0 
to N) of N groups each having two bits; 1 2a, a DA converter for converting a real-number part output from the IDFT sec- 
tion into an analog signal; 13a, an ideal filter; 14a. a multiplying section for multiplying an ideal filter output D(t) by 
cosset to perform frequency conversion. 

50 

(b-2) Arrangement of Demodulator 

When only D(t) is obtained by received frequency conversion, original information cannot be extracted without sam- 
pling 2N points. This is because D(t) consists of the real-number part of IDFT at N points. This is apparent from a sam- 
55 pling theorem as shown in FIGS. 20 A and 20 B. According to the sampling theorem, a signal in a band N • At must be 
sampled at a frequency of 1/(2 • N • Af) . Since the following equation is established: 

1/(2-N-Af) = At/2. 
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the signal must be sanopled at an interval of At/2 but an interval of At. For this reason, sampling at 2N points must be 
performed in Ts zone (N • At). 

However, when a real-number part D(t) and an imaginary-number part l(t) are obtained, an original signal can be 
extracted by sampling at N points as in the following description. 

A complex basebarxi signal obtained after received frequency conversion is expressed in Equation (5). 

[Equation 6] 

y*(0= D{0+yV{0 (5) 

The right-hand real -number part and the right-hand real-number part are obtained by deforming D(t) and l(t) by the 
transmission path and noise. When an ideal transmission path is used, the right-hand real-number part and the right- 
hand real -number part are equal to D(t) and l{t), respectively. 

'5 [Equation 7] 

N-l ^ 

y^iO — 2 [cLifoslTif^i - b^in tTtJ^t) +y(aj^m Inf^t-'r h^os 27tht) 

20 N-l ^ . 

= Z (Cjt + jb^ coslnf^t^ j{a^ -h ySjt) sin 2 jr/^ t 
Jt=o 

= Z d(Jc) e 

25 tsO 

Assume that Equation (6) is sampled by mAt (m = 0, 1 , 2,... N-1). In this case, when the following equation is estab- 
lished: 



30 
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50 



55 



[Equation 8] 



Jb=0 

= 2 dCk) e ^ 
y(m)=-^ (7) 

Equation (7) can be expressed by the following equation: 
[Equation 9] 

, N-l n ^ 



(8) 



Therefore, d(K) is obtained by Equation (8) as described follows. 
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[Equation 10] 



/ruzO 

(9) 



According ta this equation, an estimation value of the original signal d(k) is obtained. As only a reference, the relational 
JO expression between DFT and IDFT is given by the following equation: 



[Equation 11] 



N-l -a,±L 

71=0 



20 



(10) 



As described above, a complex baseband dignal obtained by frequency-converting a received signal S{t) is con- 
25 verted into a digital signal through a low-pass filter When ttie digital signal is subjected to DFT by the DFT section, the 
estimation value of the original signal d(k) is obtained. FIG. 21 is a view showing the arrangement of a main part of a 
receiver constituted while giving attention to this point. Reference numeral 15 denotes a frequency converter; 16a and 
16b, low-pass filters; 17a and 17b, A/D converters; and 18, a DFT section. 



30 (c) Frequency Conversion on Transmission Side 



An orthogonal balance modulation scheme using D(t) and l(t} is the same as a frequency conversion scheme used 
in an SSB scheme as shown in FIG. 22. Refen^ing to FIG. 22, reference numeral 11 denotes an IDFT section; 12a and 
12b, AD converters; 13a and 13, low-pass filters; and 14, a frequency converter. The frequency converter 14 is consti- 
35 tuted by multiplying sections 14a and 14b for multiplying cosco^t and sinco^t and a hybrid circuit 14c for synthesizing 
outputs from the multiplying sections to output a synthesized signal. 

An output signal S(t) from the frequency converter 14 is given by the following equation: 



[Equation 12] 

( ^if^^^ '^^tk^' ^jt^tn Infill) COS ZnJ^ (a j^in Infj^t-^ btcos Irtf^t) X ( -sin Inf^t) 

N-l 

(11) 



The signal S(t) does not include a lower sideband. Tlierefore, in this scheme, a band which is 1/2 the band in a both- 
50 sideband scheme can be obtained, arxi transmission efficiency is improved. 



(d) Frequency Conversion Scheme on Reception Side 



FIG- 23 is a view showing the arrangement of an orthogonal frequency converter using coso^t and sinoo^jt. Ref- 
55 erence numeral 1 5a denotes an adder (not actually exist] for adding a noise signal n(t) expressed by the following equa- 
tion: 

n(t) = n,t{k)cos27r(f ^ -i-f ,,)t - n Jk)sin27c{f ^ + f ,)t 
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to a received signal s(t) (see Equation (11)) to output a signal r(t); 15b, a bandpass filter; and 15c and 15d, multiplying 
sections for multiplying a bandpass filter output by cosco^t and -slnco^t . In an ideal transmission path which is free 
from amplitude attenuation, noise, and phase delay. r(t} = s(t) is established. The following description is on the 
assumption that r(t) = s(t) is established. Output signals D'(t) and r(t) from the orthogonal frequency converter are 
5 expressed by the following equations, respectively: 

[Equation 13] 

N-l 



25 



^~-D{t) (12) 

N-l 

= -f^('> (13) 

Note that, in Equations (1 2) and (13), the term of 21^ is neglected. Equations (12) and (13) coincide with the real-number 
and imaginary-number parts of Equation (3) representing a complex baseband signal on the transmission side, respec- 
tively. Therefore, as has been described above, an original signal can be extracted by performing an N-point sanrpling 
DFT arithmetic operation. 

(e) Differential Encocfing 



Although the frequency conversion scheme on reception side has been described in (d). it is very difficult to form a 
reception local frequency synchronized with a transmission local frequency a>c in the OFDM scheme. When the recep- 
35 tion local frequency includes a frequency error (asynchronous state), a denxxiulated vector is rotated as a result, and 
demodulation cannot be easily performed by the absolute phase. For this reason, information is not represented by the 
absolute phase, and the information is represented by the magnitude of phase rotation. This is called differential encod- 
ing. According to this scheme, even if there is some frequency en'or, demodulation can be performed. 



40 (e-l) Differential Encoder 



FIG. 24 is a view for explaining a differential encoder on transmission side. The logical equation of a differential 
encoder 21 is given by: 



''^ [Equation 14] 



1+7 

(14) 



50 



The differential encoder 21 is to convert data DJk) which is complex-expressed by the following equations: 



[Equation 15] 

55 

D.r(k) = A,(h)+jB,{k) d,(k) = a,(k)+jb,{h) 



into djk) by the above logical equation. 
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A differential code. 

in case of (1) [AJkj.BLCk)] = (1,1), is given by the fbllcwing equation: 

[Equation 16] 

5 

D,{k) = A,{k)+JB,{k) 



10 When this value is put into Equation (14), the following equation can be obtained: 
[Equation 17] 



20 the phase does not change. 

In case of (2) [AL.(k],B|^(k)] = (-1 ,-1), the following equation can be obtained: 

[Equa1;ion 18] 

D^k)= - 1 -j 



the phase is inverted, I.e.. the phase shifts by n. 
35 In case of (3) [AL(k),BL(k)] = (1 ,-1), the follo\wing equation is established: 

[Equation 19] 

40 

= -j • 

45 

the phase shifts clockwise by nl2, 
so In case of (4) [AL(k),BL(k)] = (-1 . 1). the following equation is established: 



55 
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[Equation 20] 
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W)= i-y" 



10 



the phase shrfts counterclockwise by id2. 

15 

(e-2) Differential Decoder 



— e 



On the basis of Equation (14}, the following equation is olstained: 
^ [Equation 21] 



I>/fc)=(l+/)- 



25 



(15) 



A differential decoder 22 is to convert data dL(k) into DL{k) according to the logical equation represented by Equation 
(15) as shown in FIG. 25. Therefore, in accordance with the cases (1) to (4) of the differential decoder, the following dif- 
ferential decoding results (1) to (4) are output. 

30 

[Equation 22] 
(1) 

in case of 



(2) 



45 



50 



55 
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In case of 



Djik)= -1 -j 



[A^k), SX*)M-1. -1) 



(3) 



In case of 



di.^{k) 




D/ik)= 



[A^k), BX«^)]=(-1. 1) 



(4) 



In case of 



df.k) 
dt.m 




= l-j 



(f) Blocks of Transmission System and Reception System 

According to the above description, a transmission system and a reception system in digital audio broadcast using 
the OFDM scheme have arrangements shown in FIGS. 26A and 26B, respectively. In a frequency converter 14 in the 
transmission system, reference numeral I4d denotes an oscillator for outputting a cos signal (costo^t) of a frequency 
fc, and reference numeral 14e denotes a phase shifter for shifting the phase of the cos signal by -90** to output -sinco^t . 
In a frequency converter 15 in the reception system, reference numeral 15e denotes an oscillator for outputting a cos 
signal (coscoct) of the frequency f^, and reference numeral 15f is a phase shifter for shifting the phase of the cos signal 
by -90** to output -sinco ^t . A cos wave and a sin wave (carrier) in the transmission system are called transmission local 
signals, and a cos wave and a sin wave in the reception system are called reception local signals. 

On the transmission side, a DAB frame is constituted by a Known phase reference symbol and M data symbols. 
Each symbol is divided into N groups each having two bits, the first and second data of each group are encoded as a 
real-number part and an imaginary-number part, respectively. The real-number part and imaginary-number part of a dif- 
ferential code are sequentially input to a Fourier transformer 1 1 , the real-number part arxJ the imaginary-number part 
output from the Fourier transformer are converted into analog signals, and the analog signals are multiplied by the cos 
and sin waves of the transmission local frequency fc, respectively. The multiplication results are synthesized to be radi- 
ated in the space. 

On the transmission side, the signal radiated in the space is received, and the received signal is multiplied by cos 
and sin waves of a reception local frequency. The multiplication results are converted into digital signals, and the digital 
signals are input to a Fourier transformer 1 8. A real-number part and an imaginary-number part output from the Fourier 
transformer are differentially decoded, and the decoded signals are sequentially output as first and second data serving 
as original data. The Fourier transformer 18 executes a Fourier transformation process on the basis of the generation 
timing of a DFT window signal. More specifically, a window signal generator (not shown) detects a null signal portion 
formed between frames to output a DFT window signal serving as a Fourier transformation execution timing of each 
symbol. The Fourier transformer 18 executes Fourier transformation on the basis of the generation timing of the DFT 
window signal. 

FIG. 27 is a timing chart for explaining a DAB frame and window signals. The start portion of the DAB frame is 
called a synchronous channel, and is constituted by a null signal portion NULL and a phase reference symbol PRS 
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(Phase Reference Symbol). A guard interval GIT of 62 ns is formed after the null signal portion NULL to reduce the 
influence of a multipath. A predeterniined number (= m) of symbols are arranged after the synchronous channel, and a 
guard interval of 62 ^s is formed before each symbol. The contents of the second-half portion of a corresponding sym- 
bol are repeatedly inserted into the guard interval. 

5 The null signal portion NULL is formed to find the start of a frame. The phase reference symbol PRS is a reference 

signal for differential decoding and receives a unique pattern (known) inherent in each frame. The PRS window is a win- 
dow arranged at a predetermined time position with reference to the detected null signal portion NULL. The DFT win- 
dow indicates the timing of a DFT arithmetic operation of each symbol. 

In the above receiver for digital audio broadcast using the OFDM scheme, the sampling frequencies of AD convert- 

10 ers 17a and 17b (FIG. 26) must be made accurately equal to the transmission rate (sampling signal frequency) fs of 
digital audio data of the transmission side. When the frequencies shift, or a phase difference occurs, the original audio 
data cannot be correctly demodulated. For this reason, the error amount of the sampling frequency must be monitored, 
and feedback control must be performed such that the error amount becomes 0. 

However, in a DAB receiver on vehicle, a multipath is generated with movement. For this reason, even if the char- 

15 acteristics of the transmission path is calculated by a signal process, the <tiaracteristics always change with the pas- 
sage of tme, and the error amount of sampling frequency cannot be easily calculated. 

SUMMARY OF THE INVENTION 

20 Therefore, it is an object of the present invention to detected an error amount of sampling frequency in the environ- 
ment wherein a multipath is generated and to perform control to make the error amount zero, thereby demodulating 
original audio data. 

According to the present invention, the above object is acehived by a receiver in digital audio broadcast using an 
OFDM scheme, comprising an RF signal demodulator for receiving a signal radiated from a DAB station into the space 

25 and multiplying the received signal by cos and sin waves of a can'ier frequency to output a baseband signal, an AD con- 
verter for converting respective multiplication results into digital data at a predetermined sampling frequency, a DFT 
window generator for outputting a DFT window signal serving as a Fourier transformation execution timing of each sym- 
bol, a Fourier transformer for performing a Fourier transformation process to the digital data for each symbol to demod- 
ulate original audio data, an impulse response arithmetic section for arithmetically operating an impulse response of a 

30 transmission path on the basis of a Fourier transformation output of the phase reference symbol, an impulse barycentric 
position arithmetic section for calculating an impulse barycentric position in the DFT window on the basis of the position 
and level of each impulse included in the impulse response, and sampling frequency control means for calculating a 
change amount of barycentric position on the basis of a deviation between a DFT window position in a previous frame 
and a DFT window position in a current frame, an impulse barycentric position in the previous frame, and an impulse 

35 barycentric position in the current frame, detecting the frequency error of the sampling pulse on the basis of the change 
amount of barycentric position, and controlling the frequency of the sampling pulse to make the error zero. 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 FIG. 1 is a model of a transmission path. 

FIG. 2 is impulse response obtained when no multipath is generated. 
FIG. 3 is a view for explaining a signal in a phase reference symbol. 
FIG. 4 is a graph showing impulse response. 

FIG. 5 is a graph for explaining impulse response depending on the presence/absence of delay 
45 FIG. 6 is impulse response obtained when a sampling frequency has an error in a state wherein no multipath is 
generated. 

FIG. 7 is an example obtained when impulse response is derived while an FFT window position is controlled in a 
state wherein a multipath is generated. 

FIG. 8 is a view for explaining FFT window position control in a multipath environment. 
50 FIG. 9 is a view for explaining a method of calculating an impulse barycentric position. 

FIG. 10 is the relationship between an FFT window position and a barycenter (when no sampling frequency error 
exists). 

FIG. 1 1 is the relationship between the FFT window position and the barycenter (when a sampling frequency error 
exists). 

55 FIG. 12 is the arrangement of a DAB receiver according to the present invention. 
FIG. 13 is a view showing the theoretical arrangement in digital audio broadcast. 
FIG. 14 is a view for explaining the function of a frequency multiplexer. 
FIG. 15 is a view for explaining a symbol. 
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FIG. 16 is a view showing the theoretical arrangement of a receiver in digital audio broadcast. 

FIG. 17 is a view for explaining D(m). 

FIG. 18 is the characteristics of an ideal fitter. 

FIG. 19 is a view showing the arrangement of a main part of a transmitter using ID FT 

FIG. 20 is a view for explaining a sampling theorem. 

FIG. 21 is the arrangement of a main part of a receiver using DFT. 

FIG. 22 is the arrangement of an orthogonal balance modulation scheme. 

FIG. 23 is the arrangement of an orthogonal frequency conversion scheme. 

FIG. 24 is a view for explaining a differential encoder. 

FIG. 25 is a view for explaining a differential decoder. 

FIG. 26 is the arrangements of a transmission system and a reception system. 
FIG. 27 is a view for explaining a DAB frame. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

(A) Principle of the Present Invention 

(a) Shift Amount of Transmission/reception Frame 

A DAB receiver specifies a shift amount between a DAB signal frame and a DAB signal frame transmitted by the 
above method. 

FIG. 1 is a model of a transmission path. Reference numeral 51 denotes a DAB signal transmitter; 52, a DAB 
receiver; and 53, a transmission path. Reference symbol x(n) denotes a transmission signal; h(n). the characteristic of 
the transmission path, and y(n). a reception signal. The respective signals have the following relationship: 

X(k).H(k) = Y(k) 

i.e., 

H(k) = Y(k)/X(k) (16) 

where X{k), H(k), and Y(k) are ot>tained by performing Fourier transformation to x(n), h(n), and y(n), respectively. 

Due to the above relationship, when x(n) and y(n) or X(k) and Y(k) are known, h(n) (impulse response) is derived, 
and a shift amount between a signal frame transmitted by the impulse response and a received signal frame can be 
specified. More specifically, a reference phase symbol (PRS), included in the start of the frame shown in FIG. 27, for 
capturing synchronization is represented by X(k). and the received PRS is represented by Y(k). thereby deriving h(n). 

(b) Relationship between Impulse Response h(n) and Shift Amount of Frames 

The relationship between the innpulse response h(n) and the shift amount of frames will be described below. When 
the sanpling frequencies of the transmitter 51 and the receiver 52 are equal to each other in the state wherein no mul- 
tipath is generated, as shown in FIG. 2, an impulse response IPL is observed at a position where n = 0 is established. 
In FIG. 2, WDW denotes a DFT window. The reason why the impulse response IPL is observed will be described below. 

In DAB, as X(k), a signal expressed by the following equation is used: 

X(k) = e^Tc^"* (-m<k<1, 1<k^m) 

0 (m does not satisfy the above conditions) This signal is illustrated in FIG. 3 (carrier exists in the range repre- 
sented by -m<k<-1 and l^k^m). Therefore, the transmission path has no delay time. H{k) is given by the following equa- 
tion: 

H{k) = Y{k)/X(k) = 1 (-m<k<-1. 1<k<m) 

0 (m does not satisfy the above conditions) In the above equation, assume that H(1) = 1 is approximately estab- 
lished. In this case, it is apparent that inverse DFT of H{k) is given by the following sine function. 

h(n) = (m/N) • [sin(n7cm/N)/(n7im/N)] (17) 
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More specifically, since m = 1536 and N = 2048 are satisfied in mode I of the DAB, the h(n) is given by h(0) = 1 where 
n = 0, and the h(n) satisfies h(n) « 1 where n ^ 0. As shown in FIG. 4, a sharp peak appears at a position where h(n] 
is maxtmum, i.e.. n (= 0). Therefore, when the transmission path has no delay time, an impulse response shown in FIG. 
5A is obtained. However, when the transmission path has a delay time corresponding to no sampling, the impulse 

5 response is expressed by h(n - no). For this reason, when n is substituted for n - no in Equation (17), the impulse 
response changes into that shown in FIG. 5B, and an impulse appears at a position which shifts by no- 

As described above, since the impulse position shifts due to signal delay, the Inrpulse position shifts even if the 
sampling frequencies of the transmitter and the receiver are not equal to each other. For example, when the sampling 
frequencies of the transmitter and the receiver are not equal to each other in the state wherein no muttipath is gener- 

10 ated, an impulse response is observed as shown in FIG. 6. 

Therefore, when no multipath is generated, the sampling frequency of the receiver is preferably controlled such that 
an impulse in the impulse response is generated at a predetermined position of the DFT window, e.g., the position 
where n = 0 is established. 

However, when a multipath is generated, as shown in FIG. 7, a plurality of impulses appear in the impulse 
15 response, and the number, positions, and levels of the impulses change with time by a change in multipath environment 
caused by movement of the receiver. In addition, the DFT window position changes to be slid forwardybackward in the 
symbol zone such that an error rate is minimized depending on the state of the multipath. FIG. 8 is a view for explaining 
a method of controlling a DFT window position. Reference symtx)l a denotes a guard interval, and reference symbol Ts 
denotes a zone of symbol B. In the guard interval a, the second-half portion of the symbol B is repeatedly inserted. The 
20 symbol position of a direct wave received by the receiver is as shown in FIG. 8A, and the symbol positions of reflective 
wave 1 to reflective wave 3 are as shown in FIGS. SB to 8D, respectively, because of the delay. As a result, a synthe- 
sized wave between the direct wave and each reflective wave is as shown in FIG. 8E. In the synthesized wave, zone a 
indicates a guard interval zone which is interfered by adjacent symbol A; zone b, a guard interval zone which is not inter- 
fered by symbol A; and zone c, the zone of symbol B. A DFT window WDW must be set at a position where each symbol 
25 is DFT-demodulated without being influenced by an adjacent symbol. The DFT window WDW is set at a position shown 
in FIG. 8E. 

As described above, in the muttipath environment, not only impulses having different positions and different levels 
exist, but also the position of the DFT window changes. Therefore, the sampling frequency of the A/D converter cannot 
be easily controlled. 

30 

(c) Outline of the Present Invention 

<1 > When sampling frequency of A/D converter has no error 

35 When a multipath is generated, an impulse indicating a direct wave cannot be easily specified because a plurality 
of impulses exist. For this reason, the barycentric position of the impulses in the DFT window are calculated, and control 
is performed on the basis of the barycentric position. As shown in FIG. 9, assume that the start position of the DFT win- 
dow is represented by 0. that the end position is represented by (N-1), that reference symtx>l N represents a sample 
number (number of data of a symbol), that reference symbol i represents a sample position in the DFT window, and that 

40 Di represents the amplitude of an impulse at a sample position i. In this case, the barycenter g of the impulses is calcu- 
lated by the following equation: 

[Equatzion 23] 

45 



(18) 

50 



Assume that the sampling frequency of the A/D converter is free from an error in a multipath generation environ- 
ment. In this case, as shown in FIG. 10A. the barycentric position of the impulses in the DFT window changes by a mul- 
55 tipath, and the DFT window position also changes as shown in FIG. 10B. For this reason, assume that the start position 
of the DFT window (PRS window) of the (L+i)th frame is represented by Xj, that the impulse barycentric position in the 
DFT window (PRS window) of the (L+i)th frame is represented by gj. and that the shift amount between the DFT win- 
dows (PRS windows) in the (L+i)th frame and the (L+i-1)th frame is represented by pj. In this case, the change amount 



N-I 

2 

1-0 
N-1 
S D, 
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Gj of impulse barycentric position between frames is given by the following equations: 

Pi = *M-^i (19) 
Gi = grgi.rPi 

Therefore, a window control amount pj between frames is given by the following equations; 

Pi = Xq - Xi 

P2 = X, - X2 



Pl = Xl-i 



In addition, the barycentric position change amount G) between frames is given by the following equations: 



^2 = g2 - gi - P2 



Gl = - gt-i - Pl 

An average value D of barycentric position change amounts between L frames is calculated by the following equation: 

D= lGj/L{i = 1 toL) (20) 

By using the average value of the barycentric position change amount, a frequency error amount AF is calculated 
by the following equation: 

AF = Dxf^3^/K (21) 

K: number of samples per one frame 

frequency of VCO (sampling frequency). 

Therefore, according to (18), (19), (20), and (21), the frequency error AF is calculated for each L frames, and the sam- 
pling frequency of the A/D converter is controlled such that the frequency error AF is 0. 

In the case wherein the sampling frequency of the A/D converter has no error as shown in FIG. 10, when the aver- 
age value D of the barycentric position change values is calculated, the following result is obtained: 

90 = 0 

gi=5. Pi=Xo-x,=n - (n-5) = 5 

Gi =91 -go-Pi =5-0-5 = 0 
g2 = -10 P2 = - Xg = (n - 5) - (n + 10) = -15 
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©2 = 92 -91 -Pi =-10-5-{-15) = 0 

/. D = 0 

Therefore, it is understood that the frequency error of the VCO for generating the sampling frequency of the A/D con- 
5 verter is 0. 

Although the above case is to control a sanipling frequency for each L frames, a sampling frequency can also be 
controlled for each frame. In this case, the frequency error amount AF is given by: 

AF = G,xf^^K (21)' 

10 

Therefore, the frequency error amount AF of each frame is arithmetically operated on the basis of (18), (19), and (21)' 
to control the sampling frequency 

(2) When sampling frequency of A/D converter has error 

15 

When the sampling frequency of an A/D converter has an error in a muHipath generation environment, as shown in 
FIG. 1 1 A, the barycentric position of impulses in a DFT window changes by a muttlpath and the frequency error. A DFT 
window position changes as shown in FIG. 1 1B. 

When the average value D of barycentric position change amounts is calculated as shown in FIG. 1 1 . the following 
20 result is obtained: 

90 = 0 

91 =7, p^ =Xo-Xi =n-(n-5) = 5 

Gi =gi -Qo-Pi =7-0-5 = 2 
25 g2 = -6 P2 - Xg = (n -5) - {n+ 10) = -15 
02=92-91 -P2 = -6-7-(-15) = 2 
.0 = 2 

Therefore, it is understood that a shift amount corresponding to two samples is generated per one frame. 
30 According to the above calculation, a plurality of impulses are generated by a multipath. Even if the piosition of the 

DFT window changes, an error component AF of the sampling frequency of the A/D converter can be extracted, and the 
sampling frequency can be controlled such that the error component becomes 0. 

(B) Embodiment of the Present Invention 

35 

FIG. 12 is a view showing the arrangement of a DAB receiver according to the present invention. 
Reference numeral 60 denotes a reception antenna; 61, an RF signal demodulator for DAB which multiplies a 
received signal by cos and sin waves of a carrier frequency to output baseband analog signals D(t} and l(t); 62, an AD 
converter for converting the baseband analog signals D(t) and l(t) into digital data D(m) and l(m) at a predetermined 

40 sampling frequency; 63, a sampling pulse generator for outputting a sampling pulse; 64, a null signal detector for detect- 
ing a null signal portion between DAB frames to output a null detection signal NDL; and 65, a DFT window generator 
for outputting a DFT window signal and a PRS window signal serving as Fourier transformation execution timings of 
symbols constituting a DAB frame. The DFT window generator 65 sets windows (DFT window and PRS window) in a 
zone which is not interfered by an adjacent symbol. 

45 Reference numeral 66 denotes a Fourier transformation/differential decoder for performing a Fourier transformation 
process to the digital data D(m) and l(m) in the PRS window and the DFT window and differentially decoding the digital 
data D(m) and l(m) to demodulate X(l^ (see Equation (16)) and audio code data (e.g., MPEG audio data); 67, an audio 
signal decoder for decoding the MPEG audio data into a PGM audio data: 68, a DA converter for converting the PCM 
audio data into analog data; 69, an amplifier; 70, a loudspeaker; 71, an operation key unit; 72, a display unit; and 73, a 

50 controller for performing tuning control or other control. 

Reference numeral 81 denotes an impulse response arithmetic section for arithmetically operating an impulse 
response by using a Fourier transformation output X(k) in a PRS window period on the basis of Equation (16); 82, an 
impulse barycenter arithmetic section for arithmetically operating an impulse barycentric position gi on the basis of 
Equation (18): and 83. a sampling error extractor for arithmetic operating an error amount AF of the sampling frequency 

55 on the basis of Equations (19) to (21) to output the error amount AF. 

The RF signal demodulator 61 for DAB multiplies a received signal by cos and sin waves of the carrier frequency 
to output baseband analog signals D(t) and l(t), and the AD converter 62 converts baseband analog signals D(t) and 
l(t) into digital data D(m) and l(m) on the basis of a sampling pulse output from the sampling pulse generator 63 to input 
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the digital data D(m) and l(m) to the null signal detector 64 and the Fourier transformation/differential decoder 66. The 
null signal detector 64 detects a null signal portion between DAB frames to output a null detection signal NDT. When 
the null detection signal NDT is input, the DFT window generator 65 moves a window to a position which is not inter- 
fered by an adjacent symbol depending on a multipath situation to input the window to the Fourier transformation/differ- 

5 ential decoder 66, and also inputs the PRS window to the innpulse response arithmetic section 81 . 

When the Fourier transformation/differential decoder 66 performs a Fourier transformation process to the digital 
data D{m) and l(m) in the PRS window and the DFT window for each symbol, and then differentially decodes the digital 
data D(m) and l(m) to demodulate X(k) and audio code data (MPEG audio data) and to output them. The audio signal 
decoder 67 decodes the MPEG audio data into PCM audio data, and the DA converter 68 converts the PCM audio data 

10 into analog data and inputs the analog data to the loudspeaker 70 through the amplifier 69, thereby outputttng DAB 
voice. 

The impulse response arithmetic section 81 arithmetically operates impulse response by using the Fourier trans- 
formation output in the PRS window period on the basis of Equation (16), arxj inputs the arithmetic result to the DFT 
window generator 65 and the impulse barycenter calculator 82. The DFT window generator 65 moves a window to a 
75 position which is not interfered by an adjacent symbol in the next DAB frame on the basis of the impulse response (mul- 
tipath). 

The impulse barycenter calculator 82 arithmetically operates an impulse barycentric position gi on the basis of 
Equation (18), and the sampling pulse extractor 83 arithmetically operates the average error amount AF of the sampling 
frequencies in L frame periods on the basis of Equations (19) to (21) to input the average error amount AF to the sam- 
20 pling pulse generator 63. The sampling pulse generator 63 controls the sampling frequencies such that the average 
error amount AF becomes 0, and inputs a sampling pulse to the AD converter 62. 

Subsequently, the above operation is repeated, the sampling frequency of the transmitter is equal to the sampling 
frequency of the receiver, and preferable audio voice can be output from the loudspeaker even in a multipath generation 
environment. 

25 In the above description, the sampling frequencies are controlled on the basis of the average en^or amourrt AF of 
the sampling frequencies in the L frame periods such that the error amount becomes 0. However, a sampling frequency 
can also be controlled for each frame. More specifically, the sampling pulse extractor 83 can also be constituted as fol- 
lows. That is, the sampling pulse extractor 83 arithmetically operates the sampling frequency error amount AF for each 
frame on the basis of Equations (1 8), (19), and (21)' and inputs the error amount to the sampling pulse generator 63 to 

30 control the sampling frequencies. 

Although the present invention has been described above with reference to the embodiment, the various changes 
of the present invention can be effected according to the spirit and scope of the invention described in the claims. The 
present invention does not exclude these changes. 

According to the present invention, the barycenter of impulse response is calculated for each signal frame of DAB, 

35 and an error conrtponent of a sampling frequency is extracted in consideration of the position of a DFT window, and the 
sampling frequency is controlled such that the error becomes 0. For this reason, the sampling frequency of a DAB 
receiver can be made equal to the sampling frequency of a transmitter even in a multipath environment, and preferable 
DAB voice can be output 

40 Claims 

1. A receiver in digital audio broadcast In which a frame is constituted by one known phase reference symtx)l and M 
data symbols, 2N digital data constituting the data symbols are divided into N groups each having two bits, the first 
data of each group and the second data of each group are sequentially input to an inverse Fourier transformer as 

45 a real-number part and an imaginary-number part, respectively, the real-number part and imaginary-number part 
output from said Inverse Fourier transformer are converted into analog signals, the analog signals are multiplied by 
a cos wave and a sin wave of a carrier frequency fc, the multiplication results are synthesized to be radiated in a 
space, signals radiated in the space are received, the received signals are multiplied by the cos wave and the sin 
wave of the carrier frequency, the multiplication results are converted into digital signals at a predetermined sam- 

50 pling frequency by an AD converter (62), the digital signals are input to a Fourier transformer (66), and a real- 
number part and an imaginary-number part output from said Fourier transformer (66) are output as the first and 
second data, 

characterized by comprising: 

55 a sampling pulse generator (63) for generating a sampling pulse to input the sampling pulse to said AD con- 

verter; 

a DFT window generator (65) for outputting a DFT window signal serving as a Fourier transformation execution 
timing of each symbol; 
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an impulse response arithmetic section (81 ) for arithmetically operating an impulse response of a transmission 
path on the basis of a Fourier transformation output of the phase reference symbol; 

an impulse barycentric position arithmetic section (82) for calculating an impulse barycentric position in the 
DPT window on the basis of the position arxJ level of each impulse included in the impulse response: and 
5 sampling frequency control means (83) for calculating a change amount of barycentric position on the basis of 

a deviation between a DFT window position in a previous frame and a DPT window position in a current frame, 
an impulse barycentric position in the previous frame, and an impulse barycentric position in the current frame, 
and controlling the frequency of the sampling pulse on the basis of the change amount of barycentric position. 

10 2. A digital audio broadcast receiver according to claim 1, characterized in that said sampling pulse generator (63) 
controls the sampling frequency such that an average error amount becomes 0. thereby inputting a sampling pulse 
to the AD converter (62). 

3. A digital audio broadcast receiver according to claim 1 . characterized in that said DFT window generator (65) sets 
IS a predetermined window at a position which is not interfered by an adjacent symbol in the next frame on the basis 

of the impulse response. 

4. A digital audio broadcast receiver according to claim 3, characterized in that the predetermined window is a DFT 
window and a PRS window. 

20 

5. A digital audio broadcast receiver according to claim 3, characterized in that the predetermined window is set at the 
position which is not interfered by the adjacent symbol when a null signal is input. 

6. A digital audio broadcast receiver according to claim 5, characterized in that the null detection signal is obtained 
25 such that a null signal portion between frames is detected by a null signal detector (64). 

7. A digital audio broadcast receiver according to daim 1 , characterized in that the impulse response indicates a mul- 
tipath situation. 

30 8. A digital audio broadcast receiver according to claim 1 , characterized in that the phase reference symbol is a PRS 
window period. 

9. A receiver in digital audio broadcast according to claim 1 , characterized in that said sampling frequency control 
means (83) controls a sanrpling frequency on the basis of an average change amount of barycentric position of L 

35 frames. 

10. A digital audio broadcast receiver according to claim 9, characterized in that said sampling frequency control 
means (83) calculates an error amount of sampling frequency by a predetermined arithmetic operation process. 

40 
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